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Abstract 

Digestion in the small intestine is the result of complex mechanical 
and biological phenomena which can be modelled at different scales. 
In a previous article, we introduced a system of ordinary differen- 
tial equations for describing the transport and degradation-absorption 
processes during the digestion. The present article sustains this sim- 
plified model by showing that it can be seen as a macroscopic ver- 
sion of more realistic models including biological phenomena at lower 
scales. In other words, our simplified model can be considered as a 
limit of more realistic ones by averaging-homogenization methods on 
biological processes representation. 
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1 Introduction 

When building a model for digestion in the small intestine, difficulties occur. The 
first one is the extreme complexity of the mechanical/biological phenomena. Trans- 
port of the bolus through the peristaltic waves, feedstuff's degradation by numerous 
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enzymatic reactions and the active/passive absorption of the nutrients by the in- 
testinal wall are known to be the key steps but they are not biologically nor fully 
understood and neither quantitavely parameterized. Modelling approaches are a 
way to integrate complex mechanisms representation of these phenomena helping 
to improve our understanding of them. Since it is almost impossible to build di- 
rect experiments for studying the digestion in the small intestine, modelling is a 
way to test in silico hypotheses that could be challenged through limited in vivo 
experiments. 

A second difficulty relies on the complex environment within the digestive tract. 
For example, the intestinal wall plays a key role in the transfer of the digested food 
in the blood and interferes in the degradation of the bolus via the brush-border 
enzymes and causes the transit of the bolus by transmitting the pulses coming 
from the peristaltic waves. 

Thirdly, digestion in the small intestine has contrasted but relevant macro- 
scopic and microscopic scales, both in space and time. To give few figures, the 
length of the small intestine in a growing pig reaches 18 meters, which is a large 
figure compared to its radius (2-3 centimeters) and even more compared to the size 
of the villi (around 1 millimeter). In the same way, the bolus stays in the small 
intestine for several hours, while the efficient peristaltic waves which ensure the 
transport of the bolus, start approximatively every 12 seconds from the pylorus. 

Because of these different scales, a model based on partial differential equa- 
tions and capturing all the interesting phenomena, would be too complicated and 
impossible to solve numerically. Therefore we have adopted in [23] a model based 
on ordinary differential equations (ode in short) : each bolus of feedstuffs coming 
from the stomach is identified as a cylinder and the odes describe the evolution of 
the position and composition of the cylinder. Since digestion could be described 
by a transport equation (or a system of such equations) with reacting terms, our 
strategy was essentially to use the Characteristics of this equation. At least nu- 
merically this type of Lagrangian method appears to be more efficient. We refer to 
[23] for details on our different models since several stages of the modelling process 
were developped in this paper. 

The aim of the present article is to provide mathematical justifications of some 
assumptions of the modelling presented in [23] . We focus on the bolus transport 
and on the effects related to absorption and enzymatic breakdown by the brush 
border enzymes, phenomena which are related to averaging/homogenization type 
processes. 

More precisely, in Section 2, we examine the effects of the pulses generated 
by the peristaltic waves. Considering that the time scale for these pulses is small 
compared to the duration of the digestion i.e. that their frequency is high, we 
rigorously establish that their effect is the same as the one of a constant driving 
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force. This result is biologically very interesting since it allows to get rid of this very 
small time scale and to do the numerical computations in a much more efficient way 
opening ways to alternative experimental approaches on digestive tract studies. 
Related and more general results on the homogenization of odes can be found in 
L. C. Piccinini[18] but we point out that our case does not fall into the scope of 
[18]- 

In Section 3, we consider the complex phenomena related to the villi and 
micro-villi : the active/passive absorption by the intestinal wall and the brush 
border enzymatic reactions. In order to study these phenomena, we introduce a 
3-d model where we focus on the boundary effects. As a consequence, the other 
phenomena are highly simplified. The lumen of the small intestine is modelled as 
a cylindrical type, periodic domain whose axis is Mei, where e\ := (1,0,0). In 
order to model the villi, this domain has an highly oscillatory boundary of order 
e _1 while its radius is of order e. In this domain, we have a system of parabolic, 
transport-diffusion equations with oscillatory coefficients for the absorbable and 
non-absorbable nutrients. The key feature is the Neumann boundary condition 
which describes the phenomena on the intestinal wall : the effects of the brush- 
border enzymes together with the active-passive absorption. 

Using homogenization method, we prove that, when e tends to 0, this problem 
converges to a 1-d system of transport-reaction equation. The key issue is to show 
how the effects of the diffusion and the degradation- absorption on the highly 
oscillatory boundary are combined in order to produce the final reaction terms. 
For the readers convenience, we provide both a formal and a rigorous proof of 
this result. The formal proof gives rather explicit formulas which can easily be 
interpreted from the biological point of view. Moreover we point out that, even if 
we are using a very simplified framework, we show that it captures the key features 
of the absorption process. 

The homogenization methods used in Section 3 are based on viscosity solutions' 
theory and in particular the "perturbed test function method" of L. C. Evans [7, 8] 
: we refer to [4] and references therein for the applications of such methods for 
problems with Neumann boundary conditions and oscillatory boundary. To our 
knowledge, it is the first time that such methods are used to obtain a convergence 
of a 3-d problem to a 1-d problem. 
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2 Transport Equation 

Peristalsis is the phenomenon in which a progressive wave of contraction or expan- 
sion (or both) propagates along a tube [16, 17, 19, 23]. The peristaltic waves are 
responsible for the fluid dynamics of the contents of the small intestine and can be 
divided into segmentation and propulsive contractions. The segmentation motion 
are responsible for mixing the bolus. Propulsive contractions are responsible for 
transporting the bolus through the small intestine. The effective peristaltic waves 
generated in the pylorus reach the bolus approximatively every 12 seconds. This is 
very small compared to the time scale of digestion phenomena which lasts several 
hours. This causes the observation of very rapid variations in the velocity of the 
bolus. 

In [23], the authors present a first simplified model of bolus transport along 
the small intestine. We use Homogenization Theory to simplify this equation to 
replace the periodically oscillating velocity by an averaged one (1). This section 
provides a rigorous mathematical justification of this transport equation. 
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Figure 1: Periodically oscillating velocity and averaged one. 



2.1 Position of the problem 

In this section, we formulate a simplified version of the transport problem. The 
small intestine is represented by the interval [0, +oo) and the position of the bolus 
at time t is given by x(t) £ [0, +oo). Roughly speaking, x(t) is the distance between 
the center of bolus and the pylorus. 

The bolus is composed of different types of nutrients, say nutrients 
1, 2, • • • , K and the quantity of nutrient i at time t is denoted by yi(t) for i = 
1, 2, • • • , K and we set y(t) := (yi (t) , j/2 (£) , • • • , Vk{P))- The variation of the dif- 
ferent yi depends on its production and degradation rate which is summarized 
through the equation 

y(t) =d(x(t),y(t)) , (1) 
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where d := (d\,d2, ■■■d n ) with di : [0, +00) x R K — > [0, +00) a Lipschitz continuous 
function. Since we are mainly interested in the transport equation in this section, 
this simple equation is written to fix ideas but also because the transport equation 
will depend on the composition of the bolus y. 

The peristaltic waves are created at the pylorus and they travel along the 
intestinal wall at a quasi-constant velocity : the average wave velocity of each 
peristaltic wave is about c ~ 7.2m/ h. These waves are periodic of period denoted 
by e <C 1 and to model them, we can say that at time t, an electric signal of size 
ip(t/e) starts from the pylorus and reaches a point x of the small intestine at time 
t + x/c. Here we assume that ip(s) = if s < and on [0, +00), ip is the restriction 
of a smooth, 1-periodic function on M. 

At time t, the bolus is at the position x(t) and is reached by the wave generated 
at time s = t — x(t)/c whose intensity is tp(s/e). we assume moreover that the 
impact of this pulse on acceleration of the bolus is given by a smooth, positive 
function g e (s, v, x, y) where, as above, s is the time when the pulse was generated, v 
is the relative velocity of a pulse with respect to the bolus velocity (v = (c—x(s))/c 
), x is the position of the bolus and y its composition. 

Indeed, according to [26] and [20] the efficiency of the peristaltic waves increases 
with the size of the bolus which is roughly speaking the sum of the y% for 1 < i < K, 
and decreases with the distance from pylorus x(t). 

The function g t is also e-periodic in s, we emphasize this fact by writing 

g e {s,v,x,y) = g(s/e,v,x,y) , 

where g(s, v, x, y) is a smooth function which is 1-periodic in s for s > 0. 

Taking into the friction inside the small intestine as in [23] through a —k(t)x(t)- 
term where k(t) > for any t, the equation for the transport of the bolus reads 

x{t) = g (e-\t - x(t)/c), 1 - x(t)/c, x{t),y(t)) - k{t)x(t) (2) 

with x(0) = and x(0) = vq where vo < c. 

Having in mind the example of a water wave in a channel, if the bolus velocity 
is the same or is close to the wave one, then obviously the peristaltic wave will 
have either no effect or at least a small effect on the bolus velocity. Translated 
in term of g, this means that g(t, 0, x,y) = and even g(t, v, x, y) = if v < 0. 
Thus there exist a smooth function j:lxKx[0, +00) xR^^ [0, +00) such that 
g(s, v, x, y) = g(s, v, x, y)v. We notice that, since we assume g to be positive, then 
g(s, v , x, y) > if v > while we have g(s, v, x, y) = if v < . 

Because of the dependence of Equation (2) on e, we denote their solutions by 
x e , y e and our aim is to study the behavior of these solutions as e tends to 0, and 
in particular the behavior of x e . 
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2.2 The Asymptotic Behavior 

We first rewrite the equation satisfied by x t ,y e . For Equation (2), we have 

x e = (1 - x*/c)g{<r\t - x e /c), 1 - x e /c, x e , y e ) - k{t)x e (3) 
while Equation (1) reads 

y e = d(x e ,y € ). (4) 

The initial conditions are 

x e (0) = 0, x'(0)=v , y e (0) = y , (5) 

where v q < c because of physiological reasons. 

In order to formulate our result, we introduce the function F(t, V, X, Y) given 

by 

F(t,V,X,Y)= [ g(s,V,X,Y)ds (6) 
J o 

and, recalling that g(s, V, X, Y) is 1-periodic for s > 0, we denote the averaged of 
F over a period by F(V, X, Y). Of course we have 

F(V, X, Y) = V, X, Y) = f 1 g(s, V, X, Y)ds, 

Jo 

and F, F are smooth functions since g is a smooth function. 
Example 2.1. In [23], the authors introduce the the function g as 

g(s, v, x, y) := ip(t - x/c) ■ v — — , 

a + bx 

for the real non-negative values co,c±,a and b. Where y = S" =1 ?/j. 

Our main result is the 

Theorem 2.1. Let (x e ,y e ) the unique solution of equations (3)-(4)-(5), then the 
sequence (x e ,y e ) e> o converges strongly in C 1 ([0, T], [0, +oo)) to (x,y) the unique 
solution of the averaged system of equations 

„ (t) = c-^(t) p{1 _ . {t)/ ^ j(t)i y{t)) _ k{t)±{t) (7) 

y(t) =d(x(t),y(t)) 
with the initial conditions 

s(0) = 0, i(0) = «o, y(0) = yo. (8) 
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The key interpretation of this result is the following : the effect of frequent 
pulses on the transport of the bolus is the same as the one obtained through an 
averaged constant signal. 

Proof of Theorem 2.1. We prove it in two steps : first we obtain various esti- 
mates showing that the sequences (x e ,y e ) €> o converge strongly in C 1 (at least 
along subsequences) and, then, in the second step, we prove that they converge to 
the unique solution of the averaged system (7) (which will imply that the whole 
sequence converges by a standard compactness argument). 

The following lemma provides a proof of convergence of x e and y e . 

Lemma 2.1. Let (x e : y e ) t> Q the unique solution of (3)-(4)-(5). Then x e ,y e are 
uniformly bounded in C 2 and therefore there exists a subsequence which is converg- 
ing strongly in C 1 and such that x e is converging in L°° weak-*. 

Proof of Lemma. We first prove that x e (t) < c. To this aim, we define the positive 
function <p{t) as follows 

/x , . x , \x e (t)-c if x e (t)-c>0 

<f>(t) = (x*(t) - c ) + = { K ' W (9) 

I otherwise, 

then multiply the both sides of equation (3) by <p(t) 

x e (x e {t)-c) + = °~ X g ^ tC ~ X , 1 - x € /c, x e , y 6 ^ (x e (t)-c) + -k(t)x e (x e (t) - c) + 

The right-hand side of this equation is negative since g(s, V, X, Y) is non negative 
if Y > 0, and k(t) > 0, therefore 

x e (t)(x e (t) - c) + < 

which is equivalent to 

l ,/ 

'(*)) < o • 



1 d ^2, 



2dt 

The function 4> 2 {t) is therefore a decreasing function. Furthermore, since vq < c, 
we have 

^ 2 (t)<^(0) = [(v -c) + } 2 = 0, 

which yields the result. 

Using the same method with (x)~ = max(— x,0), we can prove that x is a 
non-negative function. 

Gathering these information, we obtain that, for any t 

< x e {t) <ct , < x e {t) < c , 
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and therefore the sequence (x e ) t> o is uniformly bounded and equicontinuous on 
[0,T]. Using these informations and the equations for the y e , we also see that the 
y e are also uniformly bounded in C 1 (and even in C 2 ) and coming back to the x e 
equation we see also that the x € are also uniformly bounded in C 2 . 

Consequently the Arzela-Ascoli compactness criterion ensures that there exists 
a subsequence (x ej ,y €j ) which converges in C 1 . Moreover, since x e is bounded in 
L°°, we can also extract a subsequence such that x e ' converges in the L°° weak-* 
topology. □ 

We return now to the proof of Theorem 2.1. To simplify the exposure, we still 
denote by (x e ,y e ) the converging subsequence (x ej , x ej ) and we denote by (x,y) 
the limit. By inserting the Definition (6) into Equation (3), we get 

= (1 - i e /c) — (e-\t - *7c), 1 - x7c,x £ , y e ) - k{t)x\t) 

therefore, using the notation v e = 1 — x e /c and dropping most of the variables to 
simplify the expressions, we have 

x e (t) = [F(e _1 (s - x e /c), v £ , x e , y £ )] - - €±€ ^ ~ e ^f^ ~ k ^ 

and then integrate the both sides of equation over [0, t] 

r-t ft 



[ x e ds = e[ ^-(Fie-^s-xy^^^x^y^ds 
Jo Jo ds 



f* . OF dF dF ft 
- e J {ve dV +±e dX +ije dY )ds -J k±eds 

which leads to 

x e (t) -v = eF(e _1 (t - x e /c), 1 - x e /c, x e ,y e ) 

ft Q77 flZ7 «Z7" ft 



r ■ dF dF dF r 
- e J ^W + xt 8^ + yt dY )ds -J kieds 



since F(0, V, X,Y) = for any and Y G M, K . 

Now we have to let e tend to 0. First, since g is periodic, it is standard to 
prove that 

eF(e-H, V, X, Y) -»• F(V, X, Y)t , 

locally uniformly. Especially, it is easy to see that if n < t~ 1 t < n + 1, therefore 
eF(e-H, V, X, Y) -»• e[nF(V, X, Y) + 0(1)}. 
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In the same way, because of the definition of F and the regularity properties 
of g, for £ = V, X, Y we also have 

f) F f)F 
e— (e~H, V, X, Y) -> — (V, X, Y)t locally uniformly. 

As a consequence, since x e and y e are converging respectively to x and y in C 1 , 
we have also 

Cj F( e -i( s - x e ( S )/c)X( S ),x e ( S ),y £ (s)) -> - x(s)/c), 

dF 

uniformly on [0, T], where v = 1 — x/c. And the same is true, replacing F by 
dF 

and F by — . 

From these properties, it is easy to deduce that 

/"' 3F <9F /"* <9F <9F 

6 y (±e aA + y € w )ds ^J {s - x/c){± dx + yw )ds as e ^ ' 

for any t G [0,T]. 

On the other hand, v e = —x £ /c converges in the L°° weak-* topology to v and 
therefore 

rt 



r ■ dF , n, . ..dF , 



for any t G [0,T]. 

Gathering all these informations, we finally obtain 

x(t) - t; =(t - x(t)/c)F(l - x(t)/c, x(t),y(t)) 

-J o (s-x/c)(v w ds + x— +y— )ds 

- [ k(s)±(s)ds (10) 
■/ 

The right-hand side being C 1 in t, we deduce that x is a C 2 -function and 
derivating the both side of (10), we have the equation 

x(t) = (1 - x{t)/c)F{\ - x(t)/c, x{t),y{t)) - k(t)x(t). 

□ 
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3 On the Effects of Intestinal Villi 



As mentioned in the introduction, in the 1-d model of digestion presented in [23], 
we take into account the absorption effect by a simple absorption term through a 
Michaelis-Menten type nonlinearity and this can be assumed unreasonable when 
compared to the complexity of the involved phenomena. The same can be said for 
the enzymatic breakdown by the brush-border enzymes. Consequently the first 
aim of this section consists in giving some rigorous justification of these choices. 

Our effort is therefore to find an appropriate system of equations describing the 
different effects of the structure and the spatial distribution of intestinal villi on 
these key phenomena of digestion. Therefore we introduce a 3-d toy model which 
takes into account the complex geometry of the small intestine as well as all these 
boundary effects, but this implies unavoidable simplifications on the transport 
process. 

We start by a short presentation of the small intestine anatomy followed by 
introducing the three dimensional toy model of digestion. 

A large number of villi and micro-villi are present on the surface of the small 
intestine. Their role is to enlarge the digestive and absorptive area in the small 
intestine. They increase the area of the small intestine at least 500 times ([19]). 
The absorptive surface of the villi contains the brush border enzymes which are 
responsible of the final step of degradation (surfacic degradation) for some nutri- 
ents. This increase is therefore a key issue in the process of nutrients degradation 
and absorption. ([14]). 

These finger like villi are covered by epithelial cells. They consist of absorptive, 
goblet and entero-endocrine cells. The epithelial cells are produced in crypts, they 
migrate and become mature from the crypts to the tips of the villi([25]). More 
precisely, the absorption rate is also proportional to the distance of each of the 
villi from its tip. 

A microscopic observation of the small intestine surface is necessary in order to 
give realistic absorption and degradation shapes. The spatial aspect of absorption 
related to the distribution of villi and their absorption capacity is often neglected in 
modelling of digestion. In these models, absorption and degradation are modelled 
by a constant rate or a Michaelis-Menten process (see [23], Logan [15], ...). 

As shown in figure (2) the size of the period is small compared to the size of 
the unfold small intestine which is around 18 meters. We consider, for the sake 
of simplicity, that the villi are distributed periodically in the inner surface of the 
small intestine. 
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Figure 2: The different scales on intestinal anatomy relevant to our model. 



In this section, we seek a macroscopic description of digestion in the small 
intestine by taking into account all the effects of the presence of the villi in micro- 
scopic scale. Our approach is based on an asymptotic analysis, as e goes to zero. 
The absorption rate of the limit problem is said to be the homogenized absorption 
rate. 

3.1 Position of problem 

The small intestine is assumed to be an axisymmetric cylindrical tube with a 
rapidly varying cross section. In order to describe it, we first introduce an axisym- 
metric, smooth domain f2 which is confined in a cylinder of radius r > 0. More 
precisely, we assume 

{(xi, X2, x 3 ) E R 3 I x 2 = x 3 = 0} C n C {(x 1 ,x 2 ,x 3 ) 6R 3 | x\ + x\ = r 2 } . 

In addition, we assume that Cl is periodic in the xi-direction (say 1-periodic), 
namely {x\ + l,X2,x 3 ) G O if [xi,X2-,x 3 ) E f2. 

The small intestine is represented, for some < e <C 1 by the domain £l e given 

by 

Q e = enn{n>o} (ii) 
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Figure 3: A simple example of the domain Q t . The oscillations on the on the 
villi represents the microvilli. 

In this definition, the small intestine has an infinite length. However this as- 
sumption is not a real restriction, since we focus on the local absorption-degradation 
processes. Moreover, the x\ = part of the boundary corresponds to the pylorus 
and ed£l to the villi. It is worth pointing out that O e is e-periodic in the X\ di- 
rection, the parameter e characterizes the distance between the villi and thus it is 
natural to assume it to be very small. 

A simple example of Q e in cylindrical coordinates, can be the following 

Q e = {(r, 6, z) s.t I r |< e + et/>(z/e, 6), z > 0, 6 € [0, 2vr]} 

where z plays here the role of x\ and ip(z,8) is a 1-periodic function of z. 

We introduce two functions u e , v e : fl e x [0, T] — > R for describing the evolution 
of the concentration of feedstuff's in the small intestine. For x G ft e and t G 
[0,T], v e (x,t) denotes the concentration of the large feedstuffs molecules which 
are transformed into absorbable nutrients after different enzymatic reactions. The 
quantity u e (x,t) represents the concentration of produced nutrients at position x 
at each time t. 

The evolution of substrates u e and v 6 in the intestinal lumen is due to (i) their 
diffusion by Fick's law, (ii) their propagation through intestinal lumen by a given 
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velocity coming from the peristaltic waves and (Hi) the enzymatic reactions which 
transform v e to u £ both inside the intestinal lumen but also on the intestinal wall by 
the brush-border enzymes. When these reactions take place in the intestinal lumen, 
we call them volumic transformation, while we talk about surfacic transformation 
when they take place on the villi. 

The rate of the volumic reactions depends on the concentration of feedstuffs and 
also enzymes activity at time t and at x, namely ((x, t), where £ : [0, oo) x [0, T] — > 
R is a continuous, positive and bounded function. There is a limitation in the 
transformation which is described by 92 : R — > R, which is a bounded, increasing 
and Lipschitz continuous function such that <p(s) = if s < 0. These assumptions 
on £ and if are denoted by (Tl) in the sequel. 

Taking into account the three above-mentioned phenomena, the equation for 
the evolution of concentration of the non- absorbable feedstuffs molecules in the 
intestinal lumen reads 

— =Lo e Av £ -c(x 1 ,x/e,t)Dv e -C{xi,t)ip(v £ ) in n e x (0,T) (12) 
while for the absorbable nutrients, we have 

— = Xe&u e -c(x 1 ,x/e,t)Du e + ax 1 ,t)i P (v e ) in ft e x (0,T). (13) 

The first terms of the right-hand-sides of the above equations, where A de- 
notes the usual Laplacian 1 , are diffusion terms. The diffusion coefficients of large 
molecules of feedstuffs and small molecules of nutrients are denoted by ui e and Xa 
respectively. 

It is shown that, for a fixed temperature, the diffusion coefficient d is inversely 
proportional to the molecular weight, to be more precise for a spherical molecule 
we have 

in which k is Boltzmann constant, T is the intestinal temperature, \x the viscosity 
of the the intestinal liquid, p the molecule density and M the molecular mass. 
For fixed T and fi, this constant is very small because of the very small value of 
kM- x l z [14]. 

For reasons explained in section 3.3, we assume that 

io e := eu and \e ■= ex , 

, , . A , d 2 4> d 2 4> d 2 <j> 
1 li is a smooth function, A</> = — ^ + + 

0X-^ 0X^2 3 
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for some constants cu,x > 0. Since the nutrients molecules are smaller than feed- 
stuffs particles, we also have oj < x- The second terms of the right-hand sides are 
transport terms. The C 1 -function c : [0, +00) x £1 x [0, T] — > M 3 is modelling the 
velocity of substrates which comes from the peristaltic waves. The effect of the 
peristaltic waves is known to depend on the position in the small intestine and on 
time, this justifies the dependence of c(x\, X, t) on x\ and t, while the dependence 
on X = (Xi,X2, X3) takes into account the local effects at a lower scale. 

A priori the diffusion of bolus is small compared to its velocity through the 
small intestine and therefore u> e and Xe are expected to be smaller than c(x\,X, t). 

We assume the function c to satisfy the following properties : 
(CI) The function c(x±, X, t) is a Lipschitz continuous function which is 1-periodic 
in X\ and, if e\ = (1,0,0), then, for any x\ G [0, +00), X £ Q and t G [0,T], 



where P is a period in S7, say P := {X G Q ; < X\ < 1}. 

In addition to the regularity properties of c, this assumption means that the 
effect of the peristaltic waves is to move ahead the bolus in the small intestine. 
(C2) For any x\ G [0, +00), X £ Q and t G [0, T], divx(c) = where divx denotes 
the divergence operator in the X-variable only. 

This second assumption is justified by the incompressibility of the bolus at the 
microscopic level. 

(C3) For any xi G [0,+oo), X G d£l and t G [0,T], c(xi,X,t) ■ N(X) = 0, where 
N(X) denotes the outward, unit normal to dQ at X. 

This last assumption means that the velocity vector is always tangent to the 
boundary. It is worth pointing out that, if X = x/e then N(X) = n(x), therefore 
it is true both for X in Q and for x in Q e . As a consequence of this property, the 
nutrients reach the boundary only because of the diffusion effects. 

Once they reach the boundary, the large particles of feedstuffs can change 
of chemical structure because of the presence of brush-border enzymes. As we 
already mentioned above, this effect is called surfacic degradation of feedstuffs 
and the result is the production of the smaller absorbable molecules of nutrients 
u e . We assume moreover that a portion < (3 < 1 of these nutrients is absorbed 
instantaneously while the remaining part (a := 1 — /3) diffuses in the small intestine. 
The surfacic degradation is modelled by the Neumann boundary condition 



where, g is a continuous, positive and Xi-periodic function which represents the 
rate of surfacic degradation. 




= -q(xi, X)v e on dn £ x (0, T) 



(14) 
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On the boundary of the small intestine, there are two main effects for the 
nutrients u e . We already describe the first one which is a production of nutrients 
by the surfacic degradation. The second one is the active and passive absorption of 
nutrients, namely their transport across the intestinal wall to the blood circulation. 
An active process requires the expenditure of energy, while a passive process results 
from the inherent, random movement of molecules [14]. These different categories 
of absorption as well as the production of u e from v e on the boundary construct 
the boundary condition of Equation (12) 

X ~dn = _7 ?p( x i' x / e ) u£ ~ Va(xi,x/e,t)g a (u e ) + -g(xi, x/e)v e . (15) 

The functions r\ v and r\ a denote respectively the passive and active absorption 
rates. Both of them depend on the global position in the small intestine x\ and the 
local one x/e, by which we take into account the effect of the special physiology 
of the villi on the absorption rate which has been described at the beginning of 
this section. The dependence in time in the active absorption 7] a , describes the 
presence of energy at time t. The function g a governs the active absorption and 
depends on the nutrients categories. Typically, it is assumed to be the Michaelis 
Menten and therefore, it is a bounded, continuous, increasing function. 

We formulate the key assumptions on the functions rj p , rj a and g a 
(T2) The functions r) p (xi, X), rj a (xi, X, t) are bounded, continuous, positive, 1- 
periodic functions in Xi, and the function g a is a bounded, continuous, increas- 
ing function with g a (s) = if s < 0. Moreover, there exists rj > such that 
rj p (xi,X) > rj for any x\ G [0, +oo), X G ft. 

Finally, we complement the equations with the initial conditions 



which means that the small intestine is empty at time t = and by a Dirichlet 
boundary condition at x\ = 0, modelling the gastric emptying, namely 



where uq and vq are bounded continuous functions on [0, T] with ito(0) = and 
«o(0) = 0. 



i/(x,0) = 0, v e (x,0) = in ft, 



(16) 



v e (x,t) = v (t) for x x = 0, t G (0, T) 
u e (x, t) = u (t) for x x = 0, t € (0, T) , 



(17) 
(18) 
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3.2 Formal asymptotic 

In order to study the limit as e — > of the system (12)-(18), we first argue formally 
: we consider the following expansions (called ansatz) for the solutions u e and v e 

x 

u e (x, t) = u(xi,t) + eui(x±, —,t) + o(e) (19) 
v e (x, t) = v{x 1 ,t) + evi(xi, t) + o(e) (20) 

X X 

where u\{x\, —,t) and v±(xi, —,t) are 1-periodic functions in second variable. 

From now on, in order to simplify the notations, we systematically denote by X 
the fast variable x/e. On the other hand, the above system can be decoupled and 
we can first study the asymptotics of v e , namely only the initial-boundary value 
problem (12)-(14)-(16)-(17) and then use the result for studying the behavior of 
u e through (13)-(15)-(16)-(18). Since we use the same methods in both cases to 
obtain the homogenization results, we present the details only for the equation of 
nutrients u € while we only give the results for v e . 

We first plug these expressions of v e and u e into (13), and then examine the 
higher order terms in e. We find 

du d^u 1 du 

ft = x ^dx\ + " A ^i) " *)(^ e i + Dxui) + C{xi,t)ip(v) + o(l) (21) 

At this stage, we notice that the relevant choice for observing the effects of villi is 
indeed Xe = ex, for some positive constant %. With this choice, we obtain 

— = x^xui - c(x 1} X, tKg^ei + D xUl ) + C(x!,t)(p(v) + o(l). (22) 

The equation for the first corrector u\ (the "cell problem") is an equation in the 
fast variable X, i.e. for the functions X i-> u\(xi, X,t), x\,t playing the role of 
parameters. Setting 

du du 
P '' = d~x~i }t ' ei ' W^ 1 '^ (> '■= C{xi,t)(p{v), 

and substituting p and A in (22), we obtain the equation on Q, 

- X Au 1 + c(x 1 ,X,t)[p + D x u l } = X + 5 inn. (23) 

We argue in the same way for the boundary condition : plugging (19) and (20) 
into (15), we obtain 

du / ol \ 

x(T )a re 1 + D x u 1 ).n = - {r] p (x 1 ,X)u + rj a (xi,X,t)g a (u) - -g(xi,X)vj + o(l) . 

(24) 
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Using the introduced notations and recalling that N(X) = n(x), the above equa- 
tion gives 

{p+D x ui).N = -- (• np {x 1 ,X)u + r]a(xi, X, t)g a (u) - -g(x u X)v) +o(l) . (25) 

X V LO / 

Introducing the notations fi := u(x±,t) and v = v(x±, t) and 

@(xi,X,t,u,v) := rj p (x 1 ,X)u + r] a (xi,X,t)ga{u) g(x 1 ,X)v , 

the complete cell problem reads 

' -xAui + c(xi , X, t) [p + D x ui] = X + 5 inn 

1 (26) 
(p + DxU\)-N = @(x±,X, t, fi, v) on dU 

We assume that this problem has indeed a smooth solution u\ which is 1-periodic 
in X\. Recalling that is 1 periodic in the X\ direction and integrating (26) over 
a period P (remarking also that Axu\ = Ax(u\ +p ■ X)), we obtain 

(A + 5)\P\= X f p -A x («i + P ■ X)dX + J c( Xl ,X, t)[p + D xUl ])dX (27) 

where |P| denotes the Lebesgue measure of P. By using Green Formula 

-X / A x (ui + p ■ X)dX = -x {Dxui + p).Nda 
Jp JdP 

where N denotes the outward, unit normal to dP and where 

dp = (dp n on) u (dp n n) . 

We first point out that, because of the periodicity of u\ and the opposite orientation 
of the normal vector on both side of the cell 

X I (Dxui + p) ■ Nda = 0. (28) 
J(dPnn) 

On the other hand, recalling the boundary condition of (26) 

~X / (D x ui+p).Nd(j = / Q(xi,X,t, [i, v)do. 
JdPndn JdPndn 

Next we consider the c-term : by integration by parts 

/ c(xi, X, t)DxU\dX = / uic(xi, X, t).Nda — / uidivx(c)dX. 
JP JdP Jp 



Because of (C2), the last integral of the right-hand side vanishes, while, for the 
first one, we use similar argument as above : because of the periodicity properties 
of the velocity function c, the integral over dP n is (the same reasons as for 
(28)) and by (C3), it is also the case for the integral over dP n dfl. 

Gathering these informations, inserting them in (27) and recalling the defini- 
tion of 0, one gets 

(X + S)\P\= / [r) p (x 1 ,X)n+r) a (x 1 ,X,t)g a (n) g(x 1: X)u]da 

JdPndn u 

+p- J c( Xl ,X,t)dX. (29) 
In order to obtain the homogenized equation, we introduce 
c(xi,t) = J c(xi,X,t)dX, 

n P (xi) = / r )p (x 1 ,X)da (30) 

^ )= w^\L Kl elxuX)da 

where \dP D dVt\ denotes the area of the surface dP n dfl. With the notation 



a 



Q{xi,t, u, v) := Tj p (xi,X)u + r] a (xi, X, t)g a (u) + -g(xi,X)v , 

R(p) ;= lOPnm 

we get 

A = R(P)Q(x!,t, n, v) + c(xi, t) -p-S, (31) 

The one dimensional averaged equation of transport and absorption of nutrients 
is thus obtained by inserting the value of A and p in the equation (31) 

UU OUi — 

— + c(xi, t) ■ ei— = C( Xl ,t)<p(v) - R(P)@( Xl ,t, u, v) (32) 

The term R(P)@(xi, t, u, v) represents the global result of the different phenomena 
on the boundary of the small intestine : production of nutrients by surfacic degra- 
dation, active and passive absorption. The interesting feature in this term comes 
from the coefficient R(P) which measures the ratio between the large surface of 
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the villi compared to the relatively small volume of each cell. It therefore describes 
the effect of the geometry of the villi on the absorption and degradation processes. 

The term R(P)[fj p (xi)u + fj a (xi, t)g{u)\ gives an averaged value of absorption 
by intestinal wall, which takes into account the effect of villi folds as well as the 
differences between passive and active absorption. 

In the same way as for the nutrients u e , we may obtain the one dimensional 
homogenized equation for feedstuffs v 6 

dv dv 1 

— + c(x u t) • ei— = -C(x u t)<p(v) - R(P)-g( Xl )v. (33) 

In order to compare the homogenized equations (32)- (33) with the models 
presented in [23], we recall that, roughly speaking, in these models, the bolus is 
identified as a cylinder of fixed length and variable radius r, composed of a single 
feedstuff A which is transformed into an absorbable nutrient B through different 
types of enzymatic degradations. In fact, the main model is more sophisticated 
since A and B can appear under several forms (typically for A a solubilized and a 
non-solubilized form). 

Two degradation mechanisms are taken into account : a "volumic" one taking 
place inside the bolus and resulting from the action of pancreatic and gastric 
enzymes and a "surfacic" one taking place on the villi and resulting from the 
action of the brush-border enzymes. Then, once the absorbable nutrient B reaches 
the surface of bolus, hence the intestinal wall, the absorption is ensured by a 
Michaelis- Menten mechanism. Therefore, even if the above 3-d model is very 
simplified, the functions v and A have the same nature and represent the large 
particles of feedstuffs, as well as the functions u and B represent the absorbable 
nutrients. Furthermore the 3-d model described the same phenomena, at least on 
the boundary. 

Therefore, as we already mentioned it in the introduction, the above homog- 
enization process equations (32)-(33) justifies the rather simple form of the equa- 
tions presented in [23] : as long as we are just interested in "macroscopic" phe- 
nomena, it is reasonnable to describe the effects of the complex geometry of the 
villi, the different types of degradation and the absorption process by these odes 

Remark 3.1. In the above analysis, the effects of villi is summarized and measured 
by the (a priori large) R{P) -coefficient which described the consequences of their 
particular finger-like geometry. It is worth pointing out that this type of analysis 
can be used as well to understand the effects of villi in the intestinal tract but also 
the effects of micro-villi inside the villi. 
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3.3 The Rigorous Result and Proof 

We are now in position to state the rigorous result. 

Theorem 3.1. Assume that £1 is a C 2 -domain satisfying the properties described 
in Section 3.1, that (C1)-(C3) ; (T1)-(T2) holds and that uq,vq are continuous 
functions such that uq(0) = vo(0) = 0. Then the sequences (u e , v e ) e converge locally 
uniformly, as e — > 0, to the unique (viscosity) solution (u,v) of the system 



— + c(xi,t) ■ ei-Q^r =C{xi-,t)ip(v) - R(P)Q(xi,t,u,v) in Q T 
dv _ dv 1 _ 

— + c(x 1 ,t) ((x 1 ,t)ip(v)-R(P)- q(x!)v inQ T (34) 

u(0, t) = uo(t) and v(0,t) = vo(t) on dQx 
u(xi, 0) = v(xi, 0) = in [0, +00) 

where Q T = (0, +00) x (0, T) and dQ T = {xi = 0, t G (0, T)}. 

The averaged problem (34) can be seen as a simplified version of the more com- 
plicated initial-boundary value problem (12)-(18) : it is clearly easier to compute 
the solution of (34) than to take into account the complex geometry and boundary 
condition of (12)-(18). 

Proof of Theorem 3. 1 . Before providing the proof, we make some remarks about 
the existence and uniqueness of u e and v e . The system 12)-(18) is in fact decoupled 
and therefore we prove (by similar methods) the existence and uniqueness of v € 
and then of u e . 

The initial-boundary value problem (12)-(14)-(16)-(17) is a classical parabolic 
problem with Dirichlet and Neumann boundary conditions : it therefore admits 
smooth solutions. If one does not insist on proving the existence of smooth so- 
lutions, the existence and uniqueness of a viscosity solution of this problem can 
also be obtained by easier viscosity solutions arguments, using Perron's method 
(cf. [13], [6]) and comparison results ([3], [6]). Of course, the result for (13)-(15)- 
(16)-(18) follows from the same arguments. 

Applying the Maximum Principle (or a comparison result for viscosity solu- 
tions), it is easy to prove that < v e (x, t) < || vol loo in J7 e x [0, T] since and ||fo||oo 
are respectively subsolution and supersolution of (13)-(15)-(16)-(18) . In partic- 
ular, the v e 's are uniformly bounded. For the u e , the situation is unfortunately 
a little bit more complicated : since is a subsolution of (13)-(15)-(16)-(18), we 
have u e (x, t) > on Q e x [0, T] but it is not obvious at all to get an upper bound. 
For the time being, we assume that the u e 's are uniformly bounded and we will 
come back on this point at the end of the proof. 
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We provide the full convergence proof only in the case of the u e, s, the one for 
the v e being obtained by similar and even simpler argument. In this proof, because 
of the decoupling of our system, we assume that we already know that the v €, s are 
converging uniformly. 

In order to prove the convergence of u e towards u, we use the standard method 
in such problems : we combine the half-relaxed limit method [2, 6] with the Per- 
turbed Test-Function method introduced by L. C. Evans [7]. It is worth pointing 
anyway that the non-classical feature in our result and proof comes from the 3-d 
to 1-d passage to the limit and the change in the nature of the problem. 

To this end, we introduce 

u(x,t) = limsup u e (y, s) , u(x,t) = liminf u e (y,s) . 

e-*0,y-*x e-yO,y-tx 
s^t s ^ t 

We have to prove that u is a subsolution of (34) and u is a supersolution of (34); 
since the proofs for the sub and supersolution cases are similar, we only present 
the arguments for the subsolution case. 

Let <j) '■ [0, +oo) x [0, T] — > R be a smooth test-function and (x®, to) be a strict 
maximum point of u — <p. In order to prove that u is a subsolution of (34), we first 
consider the case when x\ > 0) *o > where we have to prove 

^(x?,t )+c(x?,t )-ei^-(x?,t ) < 

C(x?, to)<f>(v) - fl(P)e(sS, to, to), to)) (35) 

To apply the perturbed test-function method, we need the 

Lemma 3.1. The cell problem (26) has a Xi-periodic solution u\ if and only if 
the parameters X,p, fx, u, 5, x±, t satisfy Equation (31). Moreover this solution is 
unique up to an additive constant. 

Proof. The proof is standard and relies on the Fredholm alternative. By (C2), the 
operator C := — xA — c • D together with Neumann boundary conditions, is self- 
adjoint and the Strong Maximum Principle shows that the kernel of this operator 
only contains the constant functions. On the other hand, the computations of the 
previous subsection ensure that the right hand side of equation (26) (including the 
boundary condition) is orthogonal to the constant functions, i.e. the kernel of C, if 
and only if Equation (31) holds. Therefore this condition implies the existence of 
a solution of (26), which is C 2 by using standard elliptic regularity. This solution 
is of course unique up to an additive constant because of the structure of the 
kernel. □ 
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We pick some constant < 7 <C 1. In view of Lemma 3.1, for the choice of the 
parameters x®, to, 5 := ((x®,to)ip(v(xi,to)) 

v = v(xi,t ) , fi = u(xi, t ) - 7 , p = — (x°, t ) ex (36) 

and if we choose A given by Equation (31), there exists a smooth solution u±(X) 
of (26) associated to these parameters. 

We use this function to introduce the perturbed test-function 4> e 

X 

4> e {x,t) = 4>{xi,t) + eui(-) . 

By standard results ([2], p. 88), for e small enough, there exists a maximum 
point (x t ,t t ) of u e — near ((x°, 0, 0), to)- Moreover 

x e -> (a;?, 0,0), t e ^t as e^O (37) 

u £ {x £ ,t € ) -> u(x?,t ) as e^O (38) 

First, we prove that the maximum point (x e ,t e ) can not be on the boundary 
for e small enough. Otherwise, if (x e ,t e ) G dQ, t x (0,T), then, by the maximum 
point property on the boundary 

^(u e (x e ,t e ) - <f,(x{,t e ) - e Ul (X^)) > 

where X e = x e /e, thus 

|^V,^) - [^(xl^). ei + D xUl (X%n > 0. 

Using the smoothness of (f> and recalling that n(x t ) = N(X e ), we can write this 
inequality as 

— (x £ , f ) -]p + D XU1 (X*)} ■ N > o(l), (39) 

where, here and below, o(l) denotes a quantity which goes to as e tends to 0. 
Recalling Equation (15) permits to obtain 

~Q(xi,X e ,t e ,u e (x\,f)y(xl,t e ))-(p + D x ui(X e )j - N>o(l). (40) 

Because of (37)- (38) and the continuity properties of the functions r) p ,rj a , g and g a 
and because of the convergence of v t , Equation (40) gives 

-6 X e , t , t°), w(z°, t ))) - X (p + D x ui(X^ ■ N > o(l). 
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Now we replace u(x°, t°) by /i + 7 

(x?,X e ,*o,M + 7, *°))) -xfp + £>W* e )) -JV>o(l). 

Because of the properties of rj p , rj a , g and g a , Q(x±, X, t, u, v) is a strictly increasing 
function in u (uniformly wrt the other parameters); by using (26) together with 
(T2), we obtain 

-27 > -6 (x°, X", t , M + 7, vix^t ))) + 6 (x°, A £ , t , *°))) > o(l) (41) 

which yields the contradiction. 

Therefore the maximum point (x e , t e ) of u t — (p t is in J7 e x (0, T). Since u e is a 
solution of (13), classical properties yield to the inequality 

^(x\,f) - xAxmpTHc^X^) ^(xf,t e ) ei + £>x«i(X e 

-C(x 6 i,* e M«)<o(l). (42) 
For e small enough, using (37) and the regularity of <fi imply 

^(x?,*°) -xAxni(X e ) + c(x°, A^ )(p + ^W^)) - 5 < o(l). (43) 
Furthermore, by Equation (26), 

A = -xAui + c(x?, X e , i )(p + L>x^i) - ^ 
which yields to the inequality 

^(x?,t ) + A<o(l). 

As we already mentioned it above, this inequality is equivalent to the inequality 
(35) by just inserting the value of A from Lemma 3.1 into the above equation and 
letting e tend to 0. And the proof of this first case is complete. 

We should now consider the cases when the maximum point is achieved either 
for t = or at x\ = to complete the proof. 

For the initial condition (t = 0), a combination of the above proof and classical 
arguments shows that we have the viscosity inequality 

min{— + c(0, *) ■ ei— - C(0, t)<p(v) + R(P)0( Xl ,t, u,v),u}<0, 
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if x\ > 0, while, for x\ = 0, t > 0, one has 

min{— + c(xi,t) ■ e 1 - ((xi,t)<p(v) + R(P)Q(x 1 ,t,u,v)),u - u (t)} < 0, 

at ox\ 

and for the case x\ = 0, t = -which is a priori a particular case-, since tio(0) = 0, 
we can still use one of these inequalities which are the same. 

It is proved (cf. ([2]p.99)[6]) that, if u is a subsolution of (34), then the above 
initial condition in viscosity sense reduces in fact to a classical one 

u(x?,0)<0 on[0,+oo). 

Since c(x\,t) ■ e\ > 0, the generalized Dirichlet condition reduces also to a classical 
one (cf. [2] (cor 4.1 in p. 169)), namely 

u(0,t)<u (t) on[0,T), 

as a consequence of the fact that the characteristic are pointing outward the domain 
on the boundary. 

To conclude the proof, we invoke a (strong) comparison result for (34) : such 
result is classical and it yields u < u on [0, +oo) x [0,T], implying the desired 
convergence result. 

It remains to prove that the u e, s are indeed uniformly bounded. To this aim, we 
recall that Q is a C 2 -domain and therefore there exists a xi-periodic, C 2 -function 
d : & — > [0, oo) such that 

Dd(x) ■ N < —1 on dQ. 

Because of the particular form of fi, the function d, as well as its first and second 
derivatives, are also bounded functions. 

Now, we introduce the functions w e : Q e x [0, T] — > [0, oo) given by 

w e (x,t) = ki + k 2 t + A: 3 (e| |d| |oo - ed(x/e)) , 

for some constants k\,k2,ks > 0. We first plug these functions into the boundary 
condition (15) : using that the v e, s are bounded and that the absorption terms 
are positive, the supersolution condition is satisfied if we choose ks large enough. 
Then we consider Equation (13) : since d has bounded first and second derivatives, 
the supersolution condition is also satisfied by choosing k 2 large enough. Finally 
we choose k\ large enough to treat the boundary condition (17). 

Applying the Maximum Principle (or a comparison result for viscosity solu- 
tions) gives u e (x,t) < w e (x,t) in ft e x [0, T] and the proof is complete. □ 
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